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The possible existence of compact stars made of absolutely stable strange quark matter-referred 
to as strange stars-was pointed out by E. Witten almost a quarter of a century ago. One of the 
most amazing features of such objects concerns the possible existence of ultra-strong electric fields 
on their surfaces, which, for ordinary strange matter, is around 10 18 V/cm. If strange matter forms 
a color superconductor, as expected for such matter, the strength of the electric field may increase 
to values that exceed 10 19 V/cm. The energy density associated with such huge electric fields is 
on the same order of magnitude as the energy density of strange matter itself, which, as shown in 
this paper, alters the masses and radii of strange quark stars at the 15% and 5% level, respectively. 
Such mass increases facilitate the interpretation of massive compact stars, with masses of around 
2Mq, as strange quark stars. 

PACS numbers: 21.65.Qr; 26.60.Gj; 97.10.Cv; 97.60.Jd 



I. INTRODUCTION 

The possibility that strange quark matter (strange 
matter for short) [E H H, IU , ma de up of roughly equal 
numbers of unconfined up, down and strange quarks, may 
be the absolute ground state of the strong interaction is 
known as the strange quark matter hypothesis H H 0] • 
In the latter event objects made of strange matter- 
rang ing fro m strangelets at the low baryon number end 
[ll.[^j9ljl^.[Tll[l^[l3 J iJ.I^I| to compact stars at the high 
baryon number endgll, W, IH, EE El El & [E|-would 
be more stable than their non-strange counterparts- 
atomic nuclei and neutron stars. On theoretical scale 
arguments, strange matter is as plausible a g round state 
as the confined state of hadrons [3, i, @, Hf. Despite of 
several decades of research, both theoretical and exper- 
imental, there is no sound scientific basis on which one 
could either confirm or reject the hypothesis so that it re- 
mains a serious possibility of fundamental significance for 
physics and astrophysics [l6l |. One very striking conse- 
quence of the hypothesis is the prediction of the existence 
of a new class of compact stars-called strange (quark) 
stars d H 0, H3- The heavier members of this hypo- 
thetical family of compact stars have masses and radii 
similar to those of neutron stars. In contrast to neutron 
stars, however, strange stars would form a distinct and 



disconnected branch of compact stars and are not part of 
the continuum of equilibrium co nfig urations that include 
white dwarfs and neutron stars [ll[ EH, EI| ■ 



'Electronic address: negreiros@sciences.sdsu.edu 
T Electronic address: fweber@sciences.sdsu.edu 
* Electronic address: malheiro@it a.brl 
^Electronic address: Vladimir.Usov@wcizmann.ac.il 



If strange stars should exist in the Universe, they ought 
to be made of chemically equilibrated strange matter, 
which requires the presence of electrons inside strange 
stars. The presence of electrons plays a crucial role for 
strange stars, since they may cause the formation of an 
electric dipole layer on the surfaces of such stars lead- 
ing to huge electric fields on the order of 10 18 V/cm 
[1 Q. The situation is even more extreme if strange 
stars were made of color superconducting strange mat- 
ter [24|], which could be either in the color-flavor-locked 
(CFL) phase [H, [H [H H, HI or in the 2-flavor color 
superconducting (2SC) phase [25l [2fj|. [3fj| . In the latter 
event the electric fields on the surfaces of quark stars may 
even be on the order of 10 19 V/cm [3l| HH, depending 
on electrostatic effects, including Debye screening, and 
the surface tension of the interface between vacuum and 
quark matter [HOO!. The energy density E 2 /8ir of such 
tremendously large electric fields is on the same order of 
magnitude as the energy density of strange quark matter 
itself and, thus, should be incorporated in the energy- 
momentum tensor that is used to describe strange quark 
stars. This paper outlines how this is accomplished math- 
ematically and discusses the consequences of this extra 
contribution for the bulk properties of strange stars. 
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II. GENERAL RELATIVISTIC SELLAR 
STRUCTURE 

Our discussion is performed for spherically symmet- 
ric strange quark matter stars, whose metric is speci- 
fied by the line element (ds 2 = g vfll dx L 'dx^ , where v,p — 
0,1,2,3), 

ds 2 = e^c 2 dt 2 - e A ^dr 2 - r 2 (d8 2 + sin 2 Odcf 2 ) . (1) 

The properties of the stellar matter enter Einstein's 
field equation, G v ^ = (87r/c 4 )T l/ AI , through the energy- 
momentum tensor, T^, which for electrically uncharged 
matter is given by 



(2) 



where P is the pressure and e = pc 2 the energy density 
of strange matter, and u v its contravariant four velocity. 
The equation of state of strange matter is computed from 
the MIT bag model, 



P=(e-AB)/3. 



(3) 



for a bag constant of i? 1 / 4 = 150 MeV. This places the 
energy of strange matter at around 870 MeV, well below 
the energy per particle of the most stable atomic nucleus, 
56 Fe, as well as of infinite nuclear matter [l[. The pres- 
ence of strong electric stellar fields, as considered in this 
paper, renders the expression of the energy-momentum 
tensor Eq. ^fy significantly more complicated [H, [36[ , 

71/ = {P + pc 2 )u v u» + P8 u » 

+ ^ [f^FuI + ^6^F kl F kl ^j , (4) 

where F utl is the electromagnetic field tensor. The latter 
satisfies the covariant Maxwell equations, 



(5) 



where j v stands for the electromagnetic four-current and 
g = det{g u ^). For static stellar configurations, which are 
considered in this paper, the only non-vanishing com- 
ponent of the four-current is j . Because of symmetry 
reasons, the four-current is only a function of radial dis- 
tance, r, and all components of the electromagnetic field 
tensor vanish, with the exception of F 01 and F 10 , which 
describe the radial component of the electric field. From 
Eq. ([5]) one obtains the following expression for the elec- 
tric field, 

E(r) = F 01 (r) = l e -(*+ A )/ 2 47r f r 12 p ch e A l 2 dr\ (6) 
r J Jo 

where p c h = e*/ 2 j°(r) represents the electric charge dis- 
tribution inside the star. The electric charge within a 
sphere of radius r is given by 



Q(r) =4tt / r l2 p ch e A l 2 dr' . 
Jo 



(7) 



which can be interpreted as the relativistic version of 
Gauss' law. With the aid of Eqs. © through the 
energy-momentum tensor can be written as 
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where the electric charge is connected to the electric field 
through the relation Q 2 (r) /8?rr 4 = E 2 (r) /8tt. Substitut- 



ing Eq. 
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into Einstein's field equation leads to 
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In analogy to the electrically uncharged case, next we 
define 

GQ 2 (r) 



\r) 



1 



Gm(r) 



Equations © and dTTJ then lead to [13, 



dm 
dr 



Aur 2 



Q_dQ 

c 2 r dr 



(11) 



(12) 



where m(r) is the gravitational mass contained in a 
sphere of radius r. The first term on the right-hand-sight 
of Eq. (|12p comes from the mass-energy of the stellar mat- 
ter (quark matter in our case), while the second term on 
the right-hand-side has its origin in the mass-energy of 
the electric field carried the electrically charged quark 
star. 

The hydrostatic equilibrium equation that determines 
the global structure of electrically charged quark stars is 
obtained by requiring the conservation of mass-energy, 
0. This leads to [H 



dp 

dr 
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Q dQ 



' Airr 4 dr ' ^ 
The standard stellar structure equation of electrically 
uncharged stars, known as the Tolman-Oppenhcimcr- 
Volkoff (TOV) equation, is obtained from Eq. (US]) for 
Q — ► 0. We also note that if gravity were switched off, 
Eq. (O leads to 

V 2 p e = 4Tre 2 (n q -n e ), (14) 

with n q and n e denoting the electric charge distributions 
of quarks and electrons of strange matter, respectively. 
For a gas of free electrons one has n e — pl/3n 2 so that 
Eq. (fl4|) becomes the Poisson equation 

V 2 p e = 4Tre 2 (n q ~ pl/3n 2 ) , (15) 

which determines the electron chemical potential for 
given electric quark charge densities n q . 



III. MODELING THE ELECTRIC CHARGE 
DISTRIBUTION 

Our goal is to explore the influence of huge electric 
fields on the bulk properties of strange stars. As al- 
ready described at the beginning of this paper, the elec- 
tric charge distribution that is causing these huge electric 
fields on strange stars is located in the immediate surface 
regions of such objects. We model this distribution in 
terms of a Gaussian which is centralized at the surface 
of a strange star. This is accomplished by the following 
ansatz, 

p ch {r) = k exp(- ((r - r g )/bf) , (16) 

where b is a charge constant describing the width of the 
Gaussian, and r g is the radial distance at which the 
charge distribution is centralized. The quantity k is a 
normalization constant to be determined such that 

/+oo 
p ch (r)r 2 dr = a, (17) 
-oc 

where a is a constant proportional to the magnitude of 
the electric charge distribution. In flat space-time, a 
would be the total electric charge of the system. This 
is not the case here, however, as the electric charge Q 
depends on the metric (see Eq. ([7])). In the latter event 
the total electric charge of the system can only be com- 
puted self-consistently. Substituting Eq. (fTB|) into Eq. 
(fPT)) leads to 

8ttk = a (V^b 3 /4: + r g b 2 + ^r 2 g b/2y 1 , (18) 
which establishes a connection between k and a. 



IV. RESULTS 
A. Parameters and Boundary Conditions 

Having derived the equations that describe electrically 
charged strange stars, we now proceed to solving them 
numerically and discussing theirs solutions. We begin 
with discussing the boundary conditions of the TOV 
equation of electrically charged stars, Eq. (fT5|) . These 
are: (1) e(r = 0) = e c which specifies the energy density 
at the center of the star, (2) Q(r = 0) = which ensures 
that the electric charge is zero at the star's center, (3) 
m(r — 0) = which ensures that the gravitational mass 
is zero at the star's center, and (4) P(R) — which de- 
fines the surface of the star located at r = R. Equations 
©, ©, (O, (H3) and (HU) can then be solved for the 
mass-radius relationship of electrically charged strange 
stars shown in Fig. [TJ computed for a range of different 
values for a. The curve for which a — describes elec- 
trically uncharged strange stars. For numerical reasons 
a sample value of 0.001 km was chosen for the width b of 
the electric charge distribution on the surface of a strange 
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FIG. 1: (Color Online) Mass-radius relationship of electrically 
charged strange stars. Tick marks denote the maximum-mass 
star of each sequence, whose properties are given in Table [T] 

star. This value is many orders of magnitude greater than 
the true width of the electric charge distribution, which 
is only around ~ 10 3 fm. This particular choice, how- 
ever, has no influence on the numerical outcome, which 
can be seen by inspecting 

dQ _ r 2 a exp(-((r - r g )/bf) exp(A/2) 

dr 2(V5F6 3 /4 + r g b 2 + V^r 2 g b/2) ' [ ' 

which follows from Eqs. ([7]), (fl6|) . and (fl"8f . Equation 
din]) reveals that dQ/dr — * a(r/r g ) 2 exp(A/2)(5(r — r g ) as 
the width goes to zero, b — > 0. Hence, the dQ/dr term 
in the TOV equation (fTB")) is independent of the width of 
the electric charge distribution, provided it is sufficiently 
narrowly distributed as is the case for strange stars. 



B. Sequences of Electrically Charged Strange Stars 

Figure Q] shows the mass-radius relationship of electri- 
cally charged strange quark stars for different values of 
the electric charge constant a. The masses and radii of 
such stars change at the 10 to 15% level, depending on 
the amount of electric charge carried by the star. Table 
U lists the properties of the maximum-mass star of each 
sequence displayed in Fig. [1] It is important to note that 
the electric field strengths listed in Table |T] arc the elec- 
tric fields at the surface of the quark star (where P — > 0). 
The electron layer outside the quark star, however, nulli- 
fies the electric field, and an observer at infinity will thus 
not be able to detect it. 

The pressure profiles of electrically charged strange 
quark stars shown in Fig. [5] exhibit several peculiar fea- 
tures which are absent in compact stars made of ordinary 
neutron star matter [13, [H, H3l- First, we note that the 
interior pressure profiles of quark stars are completely 
unaffected by the electric charge layer, since the latter 
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TABLE I: Properties of electrically charged maximum-mass 
strange quark stars. The quantities R and M denote their 
radii and gravitational masses, respectively. The stars carry 
given electric charges, Q, which give rise to electric stellar 
surface fields E. 



a 


R (km) 


M (M Q ) 


2 (xl0 17 C*) E 


(10 19 V/cm) 





10.99 


2.02 








500 


11.1 


2.07 


989 


7.1 


750 


11.2 


2.15 


1486 


10.5 


1000 


11.4 


2.25 


1982 


13.5 



is located in a thin, spherical shell near the surface of 
quark stars. As a consequence, the pressure decreases 
monotonically from the stellar center toward the surface, 
as it is the case for electrically uncharged quark stars and 
ordinary neutron stars. The situation changes drastically 
in the stellar surface region, however, where the electric 
charge distribution causes a sudden and very sharp in- 
crease in pressure, as shown in the insert of Fig. [2l This 
sudden increase in pressure is a result of the appearance 
of the ultra-high electric fields and the drastic change 
in the total electric charge (dQ/dr is very high in this 
region) of the system, resulting in a significant positive 
contribution to dP/dr (sec Eq. (fTB"])). This additional 
pressure, added to the system in a region that would oth- 
erwise be the surface if the star were uncharged, creates 
a new and qualitatively different region that surrounds 
any electrically charged strange quark star. Electrically 
charged quark stars thus possess two surfaces: a "bary- 
onic matter surface" , where the surface of the uncharged 
star would be located, and an "electric surface" where 
the total pressure of the star vanishes. In between these 
surfaces is the ultra-high electric field region, which we 
will call from now on the "electrostatic layer" , since it 
only exists because of the electric field. These surfaces 
along with the electrostatic layer are shown in Fig. [3] for 
the sample quark star with a = 1000 listed in Table Q] 
The mass increase of the stars shown in Fig. Q] and Table 
U originates from the mass-energy added to these stars by 
the Coulomb field. 

Next, we discuss the electric fields of charged quark 
stars. The radial dependence of these fields is shown in 
Fig. 2] for the maximum-mass stars of Table |TJ As ex- 
pected, the electric fields exhibit a very steep increase at 
the interface between the baryonic surface and the elec- 
trostatic layer. Since the electric charge is located in 
a thin spherical layer, the electric fields quickly weaken 
with increasing radial distance, and the pressure stem- 
ming from the charge contributions drops down to zero 
at the electric surface of the star. The electric fields are 
as high as 10 19-20 V/cm. As already mentioned above, 
however, an observer at infinity would not be able to see 
these electric fields since the star is surrounded by a layer 
of electrons that nullifies the electric fields. An observer 
outside of the electric layer would thus not be able to 
detect it. 
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FIG. 2: (Color online) Pressure profiles of electrically charged 
maximum-mass strange stars shown in Fig. [1] 
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FIG. 3: (Color Online) Pressure profile for maximum-mass 
star for the a = 1000 sequence. We show the definitions 
of "Baryonic Pressure" (where the baryon pressure goes to 
zero), "Electric Pressure" where the total pressure of the star 
vanishes, and the electrostatic layer between them. 



V. CONCLUSIONS 

Almost all of the ambient conditions that characterize 
compact stars tend to be extreme. This would specifi- 
cally be the case for the electric fields carried by hypo- 
thetical compact stars made of absolutely stable strange 
quark matter (strange quark stars), which could be as 
high as 10 19 to 10 20 V/cm. In this paper, we perform 
a detailed investigation of the physical implications of 
such ultra-strong fields for the bulk properties of strange 
quark stars. The sources of these electric fields are elec- 
tric charge distributions located on the surfaces of strange 
quark stars. Depending on whether these stars are 
made of regular (i.e. non-superconducting) strange quark 
matter or color superconducting strange quark matter, 
the respective charge distributions possess very different 
physical characteristics. Color-flavor-locked (CFL) mat- 
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FIG. 4: Electric field profile inside maximum-mass stars (plot- 
ted in the relevant region). 

ter, for instance, is rigorously electrically neutral in bulk 
[281 139 ] . The situation is different for strange stars made 
of regular quark matter, which do possess electrons in 
bulk. For such stars two different electric layers would be 
formed [2j, a positively charged region of strange quark 
matter in a thin surface layer of the quark star and an 
electron layer outside of the star. The core of such stars 
has a net positive charge and the forces acting on it are 
non-zero. It is this scenario that can be described by the 
model presented in this paper. Our study shows that the 
electric charge distribution can have a significant impact 
on the structure of quark stars. For this to happen, the 
pressure (or equivalently the energy density) associated 
with the electric field needs to amount at least a few 



percent of the typical pressure P ~ 100MeV/fm 3 that 
exists inside of quark stars. Assuming that E 2 /8tt ~ 
a few xlO _2 P we thus estimate the required electric 
field strengths as E ~ xq 19 " 20 V/cm. We have shown 
that electric fields of this magnitude, generated by charge 
distributions located near the surfaces of strange quark 
stars, increase the stellar mass by up to 15% and the ra- 
dius by up to 5%, depending on the strength of the elec- 
tric field. These changes are caused by both the sudden 
increase of pressure in the surface regions of electrically 
charged quark stars, as well as by the energy density 
of the electric surface field which acts as an additional 
energy-momentum source in relativistic gravity. They fa- 
cilitate the interpretation of massive compact stars, with 
masses of around 2 M Q , as strange quark stars. 

Last but not least, another interesting feature that we 
have discovered is that the electric charge gradient term, 
dQ/dr (Eq. (TlT)|) ), which emerges in the TOV equation 
(fl~3|) of any electrically charged compact star, is inde- 
pendent of the width of the electric charge distribution, 
provided the charge distribution is sufficiently narrowly 
spread over the star. This finding is not limited to strange 
quark stars but applies to any general relativistic stel- 
lar object that carries a narrowly concentrated electric 
charge distribution. 

Acknowledgments 

This work was supported by the National Science 
Foundation under Grant PHY-55873A ARRA. Manuel 
Malheiro acknowledges financial support from CNPq and 
FAPESP. 



[1] E. Farhi and R. L. Jaffe, Phys. Rev. D 30 (1984) 2379. 
[2] C. Alcock, E. Farhi, and A. V. Olinto, Astrophys. J. 310 
(1986) 261. 

[3] C. Alcock and A. V. Olinto, Ann. Rev. Nucl. Part. Sci. 

38 (1988) 161. 
[4] J. Madsen, Lecture Notes in Physics 516 (1999) 162. 
[5] A. R. Bodmer, Phys. Rev. D 4 (1971) 1601. 
[6] E. Witten, Phys. Rev. D 30 (1984) 272. 
[7] H. Terazawa, INS-Report-338 (INS, Univ. of Tokyo, 

1979); J. Phys. Soc. Japan, 58 (1989) 3555; 58 (1989) 

4388; 59 (1990) 1199. 
[8] M. S. Berger and R. L. Jaffe, Phys. Rev. C 35 (1987) 

213. 

[9] J. Schaffher, C. Greiner, and H. Stocker, Phys. Rev. C 
46 (1992) 322. 

[10] E. P. Gilson and R. L. Jaffe, Phys. Rev. Lett. 71 (1993) 
332. 

[11] N. K. Glendenning, Ch. Kettner, and F. Weber, Astro- 
phys. J. 450 (1995) 253. 

[12] N. K. Glendenning, Ch. Kettner, and F. Weber, Phys. 
Rev. Lett. 74 (1995) 3519. 

[13] Ch. Kettner, F. Weber, M. K. Weigel, and N. K. Glen- 
denning, Phys. Rev. D 51 (1995) 1440. 

[14] Y. Zhang and R.-K. Su, Phys. Rev. C 67 (2003) 015202. 



[15] M. G. Alford and D. A. Eby, Phys. Rev. C 78 (2008) 
045802. 

[16] F. Weber, Prog. Part. Nucl. Phys. 54 (2005) 193. 

[17] N. K. Glendenning, Compact Stars, Nuclear Physics, 
Particle Phystcs, and General Relativity, 2nd ed. 
(Springer- Verlag, New York, 2000). 

[18] F. Weber, Pulsars as Astrophysical Laboratories for Nu- 
clear and Particle Physics, High Energy Physics, Cos- 
mology and Gravitation Series (IOP Publishing, Bristol, 
Great Britain, 1999). 

[19] D. Page, J. M. Lattimer, M. Prakash, and A. W. Steiner, 
Astrophys. J. Suppl. 155 (2006) 497. 

[20] A. Sedrakian, Prog. Part. Nucl. Phys. 58 (2007) 168. 

[21] B. Golf, J. Hellmers, and F. Weber, Phys. Rev. C 80 
(2009) 015804. 

[22] N. K. Glendenning, Mod. Phys. Lett. A5 (1990) 2197. 

[23] M. Malheiro, M. Fiolhais, and A. R. Taurines, J. Phys. 
G 29 (2003) 1045. 

[24] M. G. Alford, A. Schmitt, K. Rajagopal, and T. Schafer, 
Rev. Mod. Phys. 80 (2000) 1455. 

[25] K. Rajagopal and F. Wilczek, The Condensed Mat- 
ter Physics of QCD, At the Frontier of Particle 
Physics/Handbook of QCD, ed. M. Shifman, (World Sci- 
entific, 2001). 



6 



[26] M. Alford, Ann. Rev. Nucl. Part. Sci. 51 (2001) 131. 
[27] L. P. Linares, M. Malheiro, A. R. Taurines, and M. Fiol- 

hais, Braz. J. Phys. 36 (2006) 1391. 
[28] K. Rajagopal and F. Wilczek, Phys. Rev. Lett. 86 (2001) 

3492. 

[29] M. Alford and K. Rajagopal, JHEP 0206 (2002) 031. 
[30] D. Page and S. Reddy, Ann. Rev. Nucl. Part. Sci. 56 
(2006) 327. 

[31] V. V. Usov, Phys. Rev. D 70 (2004) 067301. 
[32] V. V. Usov, T. Harko, and K. S. Cheng, Astrophys. J. 
620 (2005) 915. 

[33] P. Jaikumar, S. Reddy, A. W. Steiner, Phys. Rev. Lett. 



96 (2006) 041101. 
[34] M. Alford, K. Rajagopal, S. Reddy, and A. W. Steiner, 

Phys. Rev. D 73 (2006) 114016. 
[35] S. Ray, A. L. Espmdola, M. Malheiro. J. P. S. Lemos, 

and V. T. Zanchin, Phys. Rev. D 68 (2003) 084004. 
[36] M. Malheiro, R. Picango, S. Ray, J. P. S. Lemos, and V. 

T. Zanchin, Int. J. Mod. Phys. D 13 (2004) 1375. 
[37] J. D. Bekenstein, Phys. Rev. D 4 (1971) 2185. 
[38] F. de Felice, Y. Yu, and Z. Fang, Mon. Not. R. Astron. 

Soc. 277 (1995) L17. 
[39] M. Oertel and M. Urban, Phys. Rev. D 77 (2008) 074015. 



